Mobile small polaron 
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Extending the Frohlich polaron problem to a discrete ionic lattice we study a polaronic state with a 
small radius of the wave function but a large size of the lattice distortion. We calculate the energy 
dispersion and the effective mass of the polaron with the 1/A perturbation theory and with the exact 
Monte Carlo method in the nonadiabatic and adiabatic regimes, respectively. The "small" Frohlich 
polaron is found to be lighter than the small Holstein polaron by one or more orders of magnitude. 
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A free electron interacting with the dielectric polaris- 
able continuum was studied by Pekar [0 and Frohlich ||] 
in the strong and weak coupling limit, respectively. This 
is the case of carriers interacting with optical phonons 
in ionic crystals under the condition that the size of the 
self-trapped state is large compared to the lattice con- 
stant so the lattice discreteness is irrelevant ||. The 
most sophisticated treatment of this "large" or "contin- 
uum" polaron is due to Feynman and co-workers B with 
the path-integral method, substantially extended in the 
past decade || . This treatment leads to a mass enhance- 
ment, but not to a hopping conduction or to a narrow 
polaron band. 

When the electron-phonon coupling constant A is large, 
all the states in the Brillouin zone are involved in the for- 
mation of the polaron wave function, so the polaron ra- 
dius becomes comparable with the lattice constant a and 
the continuum approximation is no longer valid. Basic 
features of the small polaron were well recognised a long 
time ago by Tjablikov ||, Yamashita and Kurosawa [Q, 
Sewell H, Holstein ||, Lang and Firsov jl0| and oth- 
ers, and are described in several review papers and text- 
books So far, analytical and numerical studies 
have been mainly confined to the Holstein model with 
a short-range electron-phonon interaction. Exact diag- 
onalization of several vibrating molecules coupled with 
one electron [|l6|Jl7|1 , variational JlJ,[li| and Monte Carlo 
calculations |20(| revealed an excellent agreement with an- 
alytical results of Holstein Q and Lang and Firsov for 
the energy of the ground state and first excited states at 
large A. Polaron mass is very large in the Holstein model, 
unless phonon frequencies are extremely high. The size of 
the region, where the small Holstein polaron is localised, 
is about the same as the size of the lattice distortion, 
each of the order of the lattice constant. Both sizes are 
almost identical also for the large Frohlich polaron, but 
much larger. 

In this Letter we study a problem of the lattice po- 
laron with a long-range Frohlich interaction pl|. This 



quasiparticle has a small (atomic) size of the electron lo- 
calization region but a large size of the lattice distortion. 
While the large Frohlich polaron is heavier than the large 
Holstein polaron, the small Frohlich polaron (SFP) turns 
out to be much lighter than the small Holstein polaron 
(SHP) with the same binding energy. We argue that 
SFPs are relevant quasiparticles in the cuprates. 

A quite general electron-phonon lattice Hamiltonian 
with one electron and the "density-displacement" type 
of interaction is given by ||,[l2|,|l5| 

H = - Wc^/Cn +};tiu)< l a{dl la d (la + 1/2) 

nn' qa 
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Here a corresponds to the different phonon modes, £ mct 
is a normal coordinate at site m, and / m a(n) is the force 
between the electron at site n and the normal coordinate 

If characteristic phonon frequencies are large compared 
to the electron kinetic energy, hio > t, (nonadiabatic 
regime) then one can apply a powerful analytic method, 
based on the Lang-Firsov canonical transformation |ic| l 
and subsequent 1/A perturbation technique. Introducing 
the phonon operators as £ mQ = I] q (umqad qa +u mqQ d q «) 

with u mqQ = ^ 1/2 (2AMw qa )^ 1 / 2 e lc ' m , N the number of 
sites, and M the ion mass, one obtains the transformed 
Hamiltonian 

H = e- S He S = - ^ cr^nC^Cn 

" £ P ^C n Cn + ^qa« Q dqa + l/2). ( 2 ) 
n qa 

Here S = EmnqafK") -1 "^"/-"^)'!.^^ - h.c, 
and E p is the familiar polaronic shift, 

e p= 2NMw 2 /"'a(0)/m'a(0)cosq(m-m / ). 

mm' qa ^ a 

(3) 
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The polaronic shift is the natural measure of the strength 
of the electron-phonon interaction. It defines the 
electron-phonon coupling constant as A = E p /zt, where 
z is the lattice coordination number. The first term in 
Eq. (^|) contains the transformed hopping integral ow, 
which depends on the phonon operators as 



exp 



/ma(n) - /mc(n') 
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At large A the hopping term in Eq. (g) can be treated 
as a perturbation. Introducing a set of N zero-order 
Bloch eigenstates ( all with the same energy — E p ) 
k,0) = TV -1 / 2 ^ n cj, exp(ik • n)|0), one readily calcu- 
lates the lowest energy levels in a crystal. Up to the 
second order in the hopping integral, the result is 



E(k) =-E p -J2 *n. 
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Here |k',n qQ ) is an excited state of the unperturbed 
Hamiltonian with one electron and at least one phonon, 
riq a is the phonon occupation number. The second term 
in Eq. (^), which is linear with respect to the bare hop- 
ping t nn ' , determines the dispersion of the polaron band 
with a band-narrowing exponent (at zero temperature) 



qa 

E [/ma(0)/m'a(0) - /ma(0)/m'«(n)] cosq(m - m'). 



The third term in Eq. (||) , quadratic in t nn > , yields a nega- 
tive almost k- independent correction of the order of 1/A 2 
to the polaron level shift. It is unrelated to the polaron 
effective mass and the polaron tunneling mobility. 

In general, there is no simple relation between the po- 
laronic shift E p and the exponent g 2 which describes 
the mass enhancement, as one can see from Eq. (@) and 
Eq. (H). We now consider the case of a single disper- 
sionless phonon mode a-> qQ = u) and the nearest-neighbor 
hopping with an amplitude t. One obtains 



^ = 2M^E/~(°). (?) 

m 

^ = = 21^ E [/-(°) - /mWAnW] ■ ( 8 ) 



The effective mass renormalisation is m* /m — e 9 , where 
m is the bare band mass and 1/m* = d 2 E(k) / 'd(hk) 2 
with k — * 0. If the interaction is local, / m (n) = Kd mn 
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FIG. 1. One-dimensional model of the small Frohlich po- 
laron on chain (x) interacting with all ions of chain (o). 



(Holstein model), then g 2 — E p /(Hlj). In general, one 
has g 2 = jEp/lhuj) with a numerical coefficient 7 = 1 — 

E m /m(0)/m(l)/E m Jm'(°). which is leSS than Unit y 

for the canonical Frohlich interaction [^2| . 

To calculate 7 explicitly we introduce one and two- 
dimensional lattice models with a long-range Coulomb 
interaction between an electron and ions (see Fig.l). The 
electron in a Wannier state on a site n of the infinite 
chain (plane) (x) interacts with the vibrations of all ions 
of another chain (plane) (o) polarised in the direction 
perpendicular to the chains. A strong coupling of car- 
riers with c-axis polarised phonons (fiu ~ 75 meV) has 
been established experimentally in Yl^CuaOg+z [p3[ . 
Because of a low c-axis conductivity and high phonon 
frequency, this coupling is not screened representing an 
example of a long-range Frohlich interaction. In this way 
our model mimics a hole on the Cu02 plane (chain x) 
coupled with the c-axis apical oxygen vibrations (chain 
o) in the cuprates. The corresponding force is given by 



/m(n) = 



1)3/2 



(9) 



Here the distance along the chains m — n| is measured 
in lattice constants a, and the inter-chain distance is 
also o=l. For this long-range interaction, one ob- 
tains E p = 1.27k 2 /(2Mw 2 ), g 2 = 0.49k 2 /{2Mhuj 3 ), and 
g 2 = 0.39E p / (hcu). The effective mass renormalisation 
is much smaller than in the Holstein model, roughly as 
ml 
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Our analytical consideration is applied if u> > t, and 
A 3> 1. To extend the results to the adiabatic case 
and to the intermediate coupling we apply a continuous- 
time path-integral Quantum Monte Carlo (QMC) algo- 
rithm, developed recently J24j] . This method is free from 
any systematic finite-size, finite-time-step and finite- 
temperature errors and allows for exact calculation of 
the ground-state energy and the effective mass of the 
lattice polaron for any electron-phonon interaction. The 
method was tested on the one-dimensional (ID) Holstein 
model which has been extensively studied by other meth- 
ods. Excellent agreement with exact diagonalization 
|l6|,^|, density-matrix renormalization group |2a] and 



2 




FIG. 2. Inverse effective polaron mass in units of 
1/m = 2ta 2 /h 2 for the one- dimensional Holstein and Frohlich 
[Eq. (0)] models. Circles: lu = 1.0 1; squares: a; = 0.5 1. 

variational jll| results was found for both the ground- 
state energy and effective mass. 

Exact polaron masses of the one-dimensional model 
defined by Eq. (||) and Fig. |[ are compared with ID 
Holstein polaron masses in Fig. |[ For both phonon fre- 
quencies Tito — 1.0 1 and 0.5 1 we found SFP to be heavier 
than SHP at small A < 1, but much lighter than SHP 
in the strong-coupling regime A > 1.5. The mass ratio 
reaches 1 order of magnitude at A = 2.75 for hco = 1.0 1 
and at A = 1.75 for Tiuj = 0.5 1. This is in accordance 
with our analytical approach in the Tiuj > t regime. Thus 
the mass ratio m* FP /m* HP is a non-monotonic function 
of A (see Fig. ||). This is a consequence of the fact that 
m FP (X) is well fitted by a single exponential function, 
exp(0.73A) for Tiuj = 1.0 1 and exp(1.40A) for hu; = 0.5 1. 
This is not so for the Holstein polaron, in which case a 
crossover between two regimes occurs at A ~ 1.5. It is 
interesting that the numerical exponents found are only 
slightly smaller than that follow from the Lang-Firsov 
transformation, exp(0.78A) and exp(1.56A) respectively. 
This shows the excellent accuracy of this transformation 
even in the intermediate region of parameters, A ~ 1 and 
TiLj/t ~ 1. Note, however, that the exact exponent devi- 
ates more and more from the Lang-Firsov approximation 
with a decreasing adiabatic ratio huj/t. This is in agree- 
ment with the exact diagonalization of a two-site model 
|jl6|| , where it was shown that the Lang-Firsov approxi- 
mation overestimates the polaron mass in the adiabatic 
regime. 

We also compared our exact QMC masses with the 
canonical weak- J2j and strong-coupling Q continuum 
polaron theory, where the bandwidth is assumed to be 
infinite. To make such a comparison meaningful we de- 
termine the Frohlich coupling constant a in such a way, 
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FIG. 3. The ratio of the effective masses of the Frohlich 
and Holstein polarons in ID. Frohlich polaron is heavier at 
small A < 1.25 but much lighter at A > 1.25. 

that the ground-state energy E of the continuum ap- 
proximation is the same as the one in our model. Then 
we calculate the continuum-case mass and compare with 
our rn* SFP (X). In the Frohlich weak-coupling regime one 
has Eq — —ahu and m* = (1 + a/6)m. This mass 
appears to be well below our m* FP for A < 1. For in- 
stance, for A = 0.5 and Hoj = t, the continuum mass 
is m* — 1.119 m while our result is m* FP = 1.422 m. 
However, in the strong-coupling regime, A > 1, the 
continuum approximation overestimates the mass. Us- 
ing Pekar's ground state energy, E = —0.1085a 2 Too, 
and mass m* — 0.021a 4 m for A = 2 and huj = t we 
find m* = 17.2 m which is much larger than our mass, 
m* F p = 4.29 m. This difference does not depend very 
much on the dimensionality of the polaron. We notice 
also that if we take into account the intermediate cou- 
pling corrections to the ground state energy of the strong- 
coupling Pekar polaron, E a = -(0.109a 2 + 2.836)?kj [||, 
a continuum polaron mass m* = 1.074 m turns out to be 
much lighter than the exact one for the same A. These 
estimates underline the crucial role of a finite bandwidth. 

To check that the light small Frohlich polaron is not 
an artifact of one dimension we calculated its mass for 
the too-dimensional (2D) version of the model (|J) and 
compared it with the 2D Holstein polaron (see Fig. Q). 
At A > 1 the mass ratio m* SFP /m* SHP (see inset) shows 
even more sharp fall than in ID. While SFP is 2.5 times 
heavier than SHP at A = 1.0, they are equal at A = 1.1, 
and SFP is 36 times lighter at A = 1.3 (at this coupling 
m SHP = 400 but rn* SHP = 11). The reason for such a 
dramatic change is the very large mass of 2D SHP. At 
the same time, the mass of SFP grows exponentially but 
smoothly, similar to the ID case. The best fit to QMC 
data is exp(1.62A + 0.19A 2 ). 

The physical reason for the small mass of SFP lies in 
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FIG. 4. Inverse effective polaron mass in units of 
1/m = 2ta 2 /h 2 for the two-dimensional Holstein and Frohlich 
[Eq. (p|)] models, hut — 0.5i. Inset: Ratio m* FP /m* HP . 



the form of electron-phonon interaction. A long-range 
interaction of the type of Eq. (g) induces a lattice distor- 
tion which undergoes less changes when the carrier hops 
to the neighboring site, than a distortion induced by a 
short-range interaction. Namely, relative changes are es- 
sential for the polaron mass. One should also emphasize 
the new type of internal structure of SFP, which is best 
understood in the extreme strong-coupling limit, A — > oo. 
In this limit the Lang-Firsov transformation is exact, and 
the polaron is localised on one site n. Hence, the size of 
its wave function is the atomic size. On the other hand, 
the lattice deformation, which is proportional to the dis- 
placement force / m (n), spreads over a large distance. Its 
amplitude falls with the distance as m — n|~ 3 in our 
model. Thus we have a new situation when the size of 
the polaron and the size of lattice deformation are very 
different. Our findings suggest to generalise the defini- 
tions of the polaron and polaron "cloud" to include this 
new possibility. 

In conclusion, we have studied the small polaron prob- 
lem with the long-range Frohlich interaction. This po- 
laron has a small (atomic) size of the wave function but a 
large size of the lattice deformation. The "small Frohlich 
polaron" . propagates in a narrow band with the effec- 
tive mass much smaller than that of the Holstein small 
polaron with the same binding energy. We argue that 
small Frohlich (bi)polarons |^] as well as large Frohlich 
(bi)polarons |2l],|28||| are relevant quasiparticles in the 
cuprates, describing holes in the C11O2 plane coupled 
with the lattice distortion by a long-range interaction. 
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